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The excitations of gluonic flux tube in a periodic lattice are examined. Monte Carlo simulations from an 
anisotropic lattice are presented and the comparison with efi'ective string models is discussed. 



1. INTRODUCTION 

It is believed that the confining regime of Yang- 
Mills theory may be described by some kind of ef- 
fective string model. The energy of the flux tube 
between static quark sources in 4d SU(3) has been 
studied in detail |1I2| : the grouping of excited 
states into bands and the corresponding energy 
gap dependence on r at large distances is sugges- 
tive of the string picture, but there is a significant 
fine structure. Recent studies in 3d Z(2), SU(2) 
and compact U(l) [H] confirmed a string forma- 
tion at large separation with less pronounced fine 
structure. In the following calculations we anal- 
yse the spectrum of the QCD periodic flux tube. 
The torelon has no fixed colour sources and so 
provides a particularly favourable theoretical en- 
vironment in which to observe the onset of string 
behaviour 



2. SIMULATION DETAILS 

On the lattice, torelons may be created using 
colour-singlet traces of ordered product of link 
matrices encircling a spatial lattice direction. Ex- 
citations are obtained by projection onto longitu- 
dinal momentum pz eigenstates and onto different 
irreps of rotations in the plane transverse to the 
flux. The energy spectrum of a QCD flux tube is 
then estimated by Monte-Carlo measurement of 
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the correlation function, 

C{t) = (^t(i)^(o)) |(vac|0|O)|2e"^°*, (1) 

where Eq is the energy of the lightest state which 
can be created by the operator (j){t). Because the 
signal decreases exponentially as t is increased, it 
is very important to obtain the asymptotic be- 
haviour of the correlator as quickly as possible. 
In order to do so we choose operators for which 
the overlap with the lightest state is as large as 
possible. We use an anisotropic lattice in which 
the temporal spacing at is much smaller than 
that in the spatial directions as to exploit the 
enhanced signal-to-noise of the correlation func- 
tion at smaller temporal separations. The center 
symmetry of the Lagrangian allows us to restrict 
our torelon operators to be positive under charge 
conjugation. The lattice symmetry of rotations 
about the winding axis z, is C4,y(E)Z{TVj for pz = 
or for Pz ^ 0. C41, contains rotations of ^ and 
the reflection in the xy plane (P-parity) and has 
flve irreducible representations: Ai,A2,Bi,B2 and 
E. Z{TZ) denotes the two-element group consist- 
ing of the identity operation and the reflection 
about the midpoint on the principal axis z (TZ- 
parity). States that are even/odd under 7^-parity 
are labeled by the subscripts g/u respectively. 

Simulations were performed on three different 
lattices of extent in the z direction of 8, 12 and 
16 tts- The lattice action was Symanzik improved 
and included a two-plaquette term, designed to 
reduce cut-off artefacts in the scalar glueball 
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Figure 1. The eight prototype lattice paths used 
in the construction of torelon operators in this 
work. 



mass. The anisotropy was set to as /at = 6 with 
flg — 0.21 fm (from rg). Linear combinations of 
the eight operators shown in Fig. ^ transform- 
ing irreducibly under the lattice symmetries were 
computed on two smeared sets of links and a vari- 
ational methods was employed on the 16 opera- 
tors to optimise overlaps. 

3. EFFECTIVE STRING THEORY 

For long QCD torelons, the two transverse os- 
cillations of the tube are the only degrees of free- 
dom forced to be low-lying Goldstone modes. 
These are the degrees of freedom of an oscillat- 
ing string. In d < 26 however, quantisation spoils 
Lorentz invariance and the string picture must 
then be regarded as an effective theory. 

An important question is whether or not suf- 
ficiently long QCD torelons can be well de- 
scribed by an effective string theory, such as the 
Nambu-Goto string in d < 26, or the Polchinski- 
Strominger proposal 0. Note that the simplest 
theories do not naturally include a number of ef- 
fects that might also be relevant, such as rigidity 
terms or interactions with the QCD bulk. 

At lowest order in an expansion about \/L, all 
these pictures give a common spectrum of states 
with energy gaps 2'kN/L. The quantum numbers 
of some of these excitations for states up to TV = 3 
are given in Table This table is restricted to 



Table 1 



Lowest string energy levels and their correspond- 
ing string and QCD states. The operator Om"^^^ 
creates an m string mode of right (left) chirality. 



Level 


State 


QCD LjJ(p,) 





|0) 


s+(o) 


1 


(af +a?T,af -<T)|0) 


n(i) 


2 


(afaf +<Ta?T)|0) 






(afaf -afaf)|0) 






(afaf ~afaf)|0) 


s-(o) 




(afaf +afaf)|0) 


sr(o) 



TRWT^fm A(2r 

((af)2 -(af )2)|0) A(2) 

afaflO) I]+(2) 

{4^ +a'^,al^ ~4^m n(2) 



~Y {al^al^ + ~af4^)\{)) A(T) 
{al^al^ -~4^af)\Q) A(l) 
(«faf +afa^t)|o) ^+(1) 

^This level includes also the states 0~ and 1* with p=l 
and states with p=3 not considered in our calculations 



those states considered in the QCD calculation. 
This spectrum can be computed for the Nambu- 
Goto action and a computation of the l/Li^ cor- 
rections in the Polchinski-Strominger scheme is 
also under investigation. The picture emerging 
from the analysis of these simple classes of ef- 
fective string theories is that the fine structure 
is predominantly dependent on the longitudinal 
momentum along the string. 

4. PRELIMINARY RESULTS 

The results from the three simulations de- 
scribed in Sec. [3 are presented in Fig. |21 The 
value L — 2 hn (determined from tq) is indicated 
by a vertical line. As can be seen in the first fig- 
ure, all states are correctly ordered according to 
their string classification, N once L > 2 fm. The 
energy gap of all levels is qualitatively in agree- 
ment with the first order string prediction 2i:N / L 
as shown in the last figure. 

The second figure shows the = 1 level is close 
to the lowest-order string prediction, although de- 
viations are resolved. The difference between the 
lowest-order string model and this excitation ap- 



pears to grow as L is increased. 

In the N = 2 states, significant structure is 
seen below 2 fm, while an extremely interesting 
pattern of degeneracies emerges at L/qs = 12 and 
16. Here, the Pz = and Pz = 2 states form two 
degenerate clusters. This is in agreement with 
indications from the effective string calculations 
that the fine structure is dependent on pz alone. 
The splitting between these clusters decreases as 
L increases. 

At iV = 3, for L > 2 fm, there appears to be 
a single degenerate set (notice from Table. ^ all 
3 states have pz — I) but a clear split between 
the Bi and B2 states is found. In the continuum 
these two states are expected to be degenerate 
and tests on possible lattice artifacts are under 
investigation. Extra operators for these states are 
being studied. 

5. CONCLUSIONS 

In this work, we have demonstrated that high- 
resolution measurements of the QCD torelon and 
its excitations can be made on the anisotropic 
lattice. This data will provide a useful forum for 
making comparisons with effective string theories 
and other models of the confining flux. Qualita- 
tive predictions of this fine structure remains a 
significant challenge. Continuum extrapolations 
are needed and simulations with larger torelons 
are also under study. 
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under Award PHY-0099450, U.S. DOE Grant 
No. DE-FG03-97ER40546, EU HPRN-CT-2000- 
00145. FM is grateful for support from Trinity 
College and Enterprise-Ireland. 
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Figure 2. Energies of torelons of different lengths. 
The solid lines are the first order string predic- 
tion, 2i:N/L. 



